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ABSTRACT 

We present detailed structure and evolution calculations for the first transiting extrasolar planets discovered by the 
space-based CoRoT mission. Comparisons between theoretical and observed radii provide information on the internal 
composition of the CoRoT objects. We distinguish three different categories of planets emerging from these discoveries 
and from previous ground-based surveys: (i) planets explained by standard planetary models including irradiation, (ii) 
abnormally bloated planets and (iii) massive objects belonging to the overlapping mass regime between planets and 
brown dwarfs. For the second category, we show that tidal heating can explain the relevant CoRoT objects, providing 
non-zero eccentricities. We stress that the usual assumption of a quick circularization of the orbit by tides, as usually 
done in transit light curve analysis, is not justified a priori, as suggested recently by Levrard et al. (2009), and that 
eccentricity analysis should be carefully redone for some observations. Finally, special attention is devoted to CoRoT- 
3b and to the identification of its very nature: giant planet or brown dwarf ? The radius determination of this object 
confirms the theoretical mass-radius predictions for gaseous bodies in the substellar regime but, given the present 
observational uncertainties, does not allow an unambiguous identification of its very nature. This opens the avenue, 
however, to an observational identification of these two distinct astrophysical populations, brown dwarfs and giant 
planets, in their overlapping mass range, as done for the case of the 8 Jupiter-mass object Hat-P-2b. According to 
the presently published error bars for the radius determination and to our present theoretical description of planet 
structure and evolution, the high mean density of this object requires a substantial metal enrichment of the interior 
and is inconsistent at about the 2-sigma limit with the expected radius of a solar-metallicity brown dwarf. Within the 
aforementioned observational and theoretical determinations, this allows a clear identification of its planetary nature, 
suggesting that planets may form up to at least 8 Jupiter masses. 
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Object 



M(Mj 



R(-Rj 



Planet 

Period a(au) | log(fi rac ) | e 



Star 

Spec. Type T eff (A') Age(Gyr) M(M Q ) [Fe/H] R(R Q ) 



CoRoT-Exo-1 b° 


1.03±0.12 


1.49±0.08 


1.50 


0.0254 


9.467 


85.1 


GOV 


5950 


0.2-4 


0.95 


-0.3 


1.11 


CoRoT-Exo-2 b b 


3.31 ±0.16 


1.465 ±0.029 


1.74 


0.0281 


9.102 


87.84 


K0V 


5625 


0.2-4 


0.97 ±0.06 




0.902 


CoRoT-Exo-3 b c 


21.66±1.0 


1.01±0.07 


4.25 


0.057 


9.277 


85.9 


F3V 


6740 


1.6-2.8 


1.37 


-0.02 


1.56 


CoRoT-Exo-4 b d 


0.72±0.08 


1 iq+o-06 

i - ly -0.O5 


9.20 


0.09 


8.468 


90 


F0V 


6190 


-i +i.o 
1 --0.3 


1.1 




1.15 


HAT-P-2 b e 


8.04±0.40 


0.98±0.04 


5.63 


0.0677 


8.962 


0.5163 90 


F8 


6290 


2.7±1.4 


1.31 


0.12 


1.48 



: Deleuil et al. (2008), d Moutou et al. (2008), ' Winn et al. 

Table 1. Transiting Exoplanets characteristics. For all the CoRoT objects, e = is assumed in the analysis of the 
light-curve and of the radial velocity data. 



2. Model Description 

The main physics inputs used in the present models have 
been described in previous works devoted to the evo- 
lution of b rown dwarfs and giant planets (jBaraffe et al.l 
120031 120081 ) . Evolutionary models are based on a consistent 
treatment between non-grey irradiated atmospheric struc- 
tures and inner structures. The interior is composed of a 
gaseous envelope of solar mixture described by the Sau mon- 
Chabrier equation of state (EOS. ISaumon et al.lll995T ). The 
inner structure models also take into account enrichment of 
heavy material, with the appropriate EOSs ([Baraffe et al.l 
2008). For sake of simplicity, the effect of heavy element 
enrichment is mimicked by the presence of a water core 
since water, with silicates, is the most abund ant species 
in the protoplanetary disk. iBaraffe et alj (|2008l ) show that 
for enrichment less than 15% by mass, as considered in 
the present study, the differences in radius between a pure 
water or a pure silicate core composition, as well as be- 
tween a heavy material centrally condensed or distributed 
throughout the envelope, are small (< 2%). We use the 
ANEOS EOS for water ([Thompson &; Lausonl [1972] ). A 
detailed analysis of the effect of the different EO S avai l- 
able in the literature can be found in IBaraffe etaTl (|2008f) . 
Co mparison with mode ls from other groups, in particular 
bv lFortnev et all (|2007l ). can also be found in the same pa- 
per. 

Transiting planets are very close to their host star 
[a < 0.1 AU) and the effect of stellar irradiation is an im- 
portant process to be accounted for in their modelling. For 
a 1 Gyr-old Jupiter-like planet orbiting a solar like star at 
0.05 AU, the incident stellar flux is more than 10 000 times 
larger than the thermal flux of the planet itself. In this 
case, stellar irradiation strongly affects the planet a tmo- 



shows several evolutionary tracks for CoRoT-4b in a radius 
versus time diagram. 



spheric structure to deep levels ([Barman et al. 2001) an d 



thus its evolution (Baraffe ct al. 20031. iBurrows et al.| [2003). 
The presence of an irradiated atmosphere acts as a buffer 
which slows down the release of the internal gravo-thermal 
energy, and thus th e cooling and contraction of the planet 
([Guillot et al.lll996fl . We use a grid of irradiated atmospher e 
models based on the calculations of iBarman et a l. (2001), 
computed for different levels of stellar irradiation relevant 
to the present study. 



3. Standard planets 

In this section we focus on planets whose radius can be ex- 
plained by the " standard" models discussed in iJ21 CoRoT- 
Exo-4b belongs to this population. The planetary and stel- 
lar properties of this system are given in Table [TJ Fig. [1] 




l 2 

Age (Gyr) 

Fig. 1. CoRoT-Exo-4b: evolution of the radius as a function 
of age. Solid line: standard cooling sequence of an isolated 
0.72Mj gaseous sphere with solar composition. Long-dashed 
line: irradiated case. Dash-dotted line: Irradiated case with 
a 10 M ffl water core. Red box: observational la error bar. 



In order to illustrate the importance of irradiation ef- 
fects in the modelling of transiting planets, Fig. [1] displays 
the evolution of an isolated gaseous sphere with solar com- 
position (solid line). In this case, the predicted radius is too 
low compared to the observed one. When taking thermal 
effects of the stellar impinging flux on the planet's atmo- 
sphere into account (long-dashed line), the model yields a 
radius consistent with observations. For this class of ob- 
jects, one can infer an upper limit for the content of heavy 
material by determining the maximum mass of heavy ele- 
ments consistent with the lower observational error bar. In 
the case of CoRoT-Exo-4b, models with less than 10 M ffi of 
water remain within the observational error box, as illus- 
trated by the dash-dotted curve in Fig. [1] This corresponds 
to a total mass fraction of heavy elements Z < 5%. As men- 
tioned in Sj2[ adopting rock as the main heavy element or a 
mixture of water and rock will barely change this value. For 
this type of transit planet population, the bulk composition 
can thus be constrained with reasonable accuracy. 

One may argue that since a physical process is miss- 
ing in the modelling of some close-in planets - as discussed 
in the section below - this upper limit for heavy material 
enrichment could be underestimated, if such a process oc- 
curs in all transiting planets. An additional heat source, 
for instance, will yield a larger planet's inflation at a given 
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age, allowing for a larger maximum amount of heavy mate- 
rial consistent with the observational error box. The exact 
nature of such additional heating mechanisms (discussed 
in 52} , or of other physical mechanisms invoked to explain 
the abnormally large radii, however, is still debated and 
may depend on the planet/star system parameters. Under 
such conditions, the wisest (and simplest) assumption is 
not to invoke missing physics when it is not required. One 
should bear in mind, however, that, given the remaining 
limitations in our present understanding of planet struc- 
ture and evolution, these heavy element enrichment deter- 
minations retain some degree of uncertainty. This maximum 
mass should rather be seen as the maximum enrichment 
consistent with the observed la error bar, according to the 
present theoretical models. 

4. Abnormally inflated planets 

A significant fra ction of transiting plan ets shows abnor- 
mally large radii (|Mandushev et al.ll2007l ) compare d to pre- 



dictions of standard irradiated models (see e.g. Guillot 
120081 IChabrier et all l2009t iBaraffe et all 120091 and refer- 
ences therein). So far, two objects in the CoRoT sample be- 
long to this category : CoRoT-Exo-lb and -2b. Comparison 
between models and observations for these two planets are 
shown in Fig. [5] and [3J respectively. The effect of irradia- 
tion alone does not allow to reproduce the observed radii 
of these strongly inflated objects. These observations thus 
confirm the existence of a missing mechanism in the mod- 
els, which slows down the cooling and the contraction of 
the planet. 



1.6 - 



& 1.4 
Pi 

1.2 - 



1.0 



1 2 3 4 5 

Age (Gyr) 

Fig. 2. CoRoT-Exo-lb. Solid line: isolated solar composi- 
tion giant planet. Long dashed line: irradiated case. Dash- 
dotted line: effect of tidal heating computed with e = 0.02. 
Short dashed line: solar composition case with 1% of the 
stellar flux transported deep into the planet interior . 

Several possibilities have been suggested to explain this 
puzzling property. Tidal heating due to circularization o f 
the orbit, as originally suggested (|Bodenheimer et aT1l2001[ ) 
and then rejected on the basis of a too short characteris- 
tic timescale compared with the age of the systems, might 
provide or at least participate to the lacking mechanism in 
some cases, provided tidal effects in the planet and in the 
star are properly taken into account (Jackso n et al. 2008) 
(jLevrard et all 120091 : llbgui k. Burrows! [20091. In a recent 
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Age (Gyr) 

Fig. 3. CoRoT-Exo-2b. Same legend as in Fig. 2. 



paper, however. iLevrard et al.l (|2009h show that most tran- 
siting systems are not in a state of tidal equilibrium and 
thus that tidal circularization timescale estimates based on 
equilibrium ti d es are not correct. A major consequence of 
ILevrard et alj (|2009l ) calculations is that assuming zero- 
eccentricity, when not directly determined, in transiting 
light curve analysis, is not necessarily correct and such anal- 
ysis should be redone carefully. Tidal heating might thus in 
some cases provide an extra source of energy during a sig- 
nificant fraction of the planet's evolution. In the present 
mod e ls, we account for this source of energy, according to 
iHutl (|1981| ). Note that, since zero eccentricity is assumed 
for the analysis of all the CoRoT transit light-curves, we 
assume a small but finite value for e in order to have a 
quantitative estimate of the effect of tidal heating on the 
bloated CoRoT planets. In the case of CoRoT-Exo-lb, an 
eccentricity e = 0.02 provides enough tidal heating to repro- 
duce the observed radius within the error bars (dash-dotted 
curve in Fig. [5]). For CoRoT-Exo-2b, an eccentricity of only 
a few percents does not provide enough tidal heating and 
larger values are required to reproduce the radius. We find 
that a value e = 0.15 yields a radius in agreement with 

ob servations. 

iShowman fc Guillo t (2002) suggested another heating 
mechanism in the deep interior of planets, originating from 
the strong winds generated at the planet's surface. Their 
numerical simulations of atmospheric circulation produce a 
downward kinetic energy flux of about 1% of the absorbed 
stellar incident flux, which dissipates in the interior and 
slows down the planet's contraction (|Guillot fc Showman 
|2002[ ) . Although the validity of this scenario is still debated, 
with various simulatio ns producing different results (see e.g. 
Sho wman et al.l 12008). it is worth exploring this issue. To 
test this effect, we have included an extra source of energy 
corresponding to 1% of the stellar impinging flux contribu- 
tion in our evolutionary models, as done in IChabrier et all 
(200|). In the case of CoRoT-Exo-lb (dashed curve in Fig. 
[2]), the resulting theoretical radius is too large, by ~ 8%, 
compared to the observed value, while for CoRoT-Exo-2b, 
the same relative contribution yields too small of a radius. 
Note that for CoRoT-Exo-4b, analysed in |j3j models in- 
cluding 1% of the stellar impinging flux as an extra source of 
energy also overestimate by 4% the radius compared to the 
observed value, except if the heavy element content of the 
planet is significantly larger than the inferred ^5% mass 
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frac tion. Enhanced atm ospheric opacity in transiting plan- 
ets (|Burrows et al.l [20071 ). although not excluded, remains 
so far too much of an ad-hoc suggestion to be examined in 
detail (see lBaraffe et al]l2009l for a discussion). At any rate, 
the presently detected most inflated transiting planets, like 
Tres-4b and WASP-12, can n ot be explained even with 
such enhanced-opacity models (|Guillotl [2008: Ba raffe et alJ 
12009ft . 

An other mechanism, based on (inefficient) layered or 
oscillatory co nvection in some pla n et int eriors, has been 
suggested by IChabrier fc Baraffel (|2007l ) and has been 
shown to possibly explain the abnormally large radii. If fu- 
ture follow-up observations of CoRoT-Exo-lb and CoRoT- 
Exo-2b confirm eccentricity values e < 0.01, such a mecha- 
nism will have to be considered with serious attention. 



5. Massive substellar objects 

The discovery of "super" Jupiters, with masses > 10Mj, in 
close orbit to a central star, raises questions about their 
nature: planet or brown dwarf ? CoRoT-Exo-3b (see Table 
[1]) is a perfect example of such an ambiguity. Studies of low 
mass stars and brown dwarfs in young cl usters suggest a 
conti nuous mass function down to ~ 6 Mj (jCaballero et alJ 
|2007[ ). indicating that the same formation process respon- 
sible for star formation can produce objects down to a 
few Jup iter masses. Indeed, analytical theories of star for- 
mation (|Padoan fc Nordlundl[2004lHennebelle fc Chabrierl 
(2008;:)) show that gravoturbulent fragmentation of molec- 
ular clouds produces, with the same underlying processes, 
stars and brown dwarfs down to a few Jupiter-masses in 
numbers comparable to the observationally determined dis- 
tribution. Brown dwarfs and planets thus overlap in mass, 
stressing the need for identi?cation criteria enabling the 
distinction between these two types of astrophysical bod- 
ies. The presence of strongly non-solar atmospheric abun- 
dances, as observed in the atmosphere of the giant planets 
of our Solar System, may provide signatures of a planetary 
formation process in a proto-planetary disk. Such a signa- 
ture, however, is difficul t both to observe and to character- 
ize at the present time (|Chabrier et alj |2007) and may not 
apply to irradiated planets, with radiatively stable outer 
layers. A more robust signature of the planet formation 
process, as expected from the core accretion model, is the 
presence of a significant amount of heavy material in the 
interior. Observed radii signi?cantly smaller than predicted 
for solar or nearly-solar metallicity objects reveal the pres- 
ence of such a signi?cant average amount of heavy material; 
a major argument in favor of the core-accretion planet for- 
mation process. On the opposite, if a physical mechanism 
is missing in current planet cooling models, as discussed in 
21 observed radii larger than predicted do not necessarily 
imply an absence or a small amount of heavy material. For 
such cases, the nature of the object remains ambiguous, if 
only based on the knowledge if its mean density. 

In this section, we focus on the case of CoRoT-Exo-3b 
and examine whether its radius determination enables us 
to identify its very nature. As shown in Fig. 01 the observed 
radius of CoRoT-Exo-3b can be matched by the model of an 
irradiated brown dwarf of 21.6 Mj with solar composition 
(long-dashed line). This is by itself an encouraging confir- 
mation of the theoretical prediction of the age-mass-radius 
relationship in the brown dwarf regime (IChabrier fc Baraffel 



12 3 4 




Age (Gyr) 

Fig. 4. CoRoT-Exo-3b. Solid line: standard cooling se- 
quence of an isolated 21.66Mj brown dwarf with solar 
composition. Long dashed line: irradiated case. Dash-dotted 
line: Irradiated case with a 272 M ffi core of water. Red box: 
observational la error bar. 

119971: iBaraffe et al.lll998t ). Note that, given the small orbit- 
ing distance, the effects of irradiation are not negligible, 
even for such massive objects. Accounting for irradiation 
on the atmospheric profile, and thus on the object's cooling 
history, is thus mandatory to provide consistent compari- 
son between models and observations, when the radius is 
determined at this level of accuracy (~ 7%). The present 
radius error bars, however, are still too large to infer or 
exclude the presence of a significant amount of heavy ma- 
terial in the interior of this object. As done in Sj3j a maxi- 
mum amount of heavy material can be determined, for the 
minimum theoretical radius allowed by the error box. We 
find an upper mass limit for the core of about 8OOM0, 
i.e. a global maximum mass fraction Z < 12%. We stress 
that this corresponds to the maximum enrichment com- 
patible with the actual error bars. The possibility to have 
such a large amount of heavy material must be examined 
in the context of our current understanding of planet for- 
mation, within the framework of the core-accre t ion m odel. 
Following our previous analysis (Bar affe et al.l (|2008| ) and 
references therein), we can estimate the maximum amount 
of heavy material available in the proto-planetary disk for 
planet formation. According to cu rrent models of pl anet 
formation which include migration (Alib ert et al.|[2005l ). up 
to 77 ~ 30% of heavy elements contained in the protoplan- 
etary disk can be incorpora ted into forming giant planets 
(jMordasini et al.ll2008l [2009). The maximum mass of avail- 
able heavy material that can be accreted to form planets is 
thus: 

M z « 77 ■ Z • (/ ■ AQ (1) 

where /• M* is the maximum mass for a stable disk (< 
0.1M*) and Z is the metal mass fraction of the star. For 
CoRoT-Exo-3b, which is orbiting a 1.37 M© F star with 
near solar metallicity, at most Mz ~ 270M ffi of heavy 
material can thus be accreted to form the planet. This 
(admitedly crude) upper limit derived from current planet 
core accretion formation models yields a planet contrac- 
tion consistent with today' s observations, as see n in Fig. 0J 
Note that, as discussed in IBaraffe et all 1)20081 ) for HAT- 
P-2b, this heavy material does not need to be accreted 
into one single object, as very massive planets, in particular 
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short-period ones, may result from smaller planet collisions. 
Therefore, given the present uncertainties in the radius de- 
termination, neither the brown dwarf nor the planet possi- 
bility can be assessed or excluded for CoRoT-Exo-3b, whose 
nature remains ambiguous. A comparison between the pre- 
dicted radius of a (irradiated) solar-metallicity brown dwarf 
(dashed line) and of a planet with the aforedetermined mas- 
sive core, which represents only a ~ 4% metal enrichment 
(dash-dotted line) in Fig. [4] shows that a radius accuracy 
< 3% is required to resolve the ambiguity, according to the 
present models. In any event, this demonstrates the promis- 
ing powerful diagnostic provided by mass-radius determina- 
tions to distinguish massive planets from low-mass brown 
dwarfs, providing adequate observational accuracy. 

Among the few known massive planetary-mass objects, 
there is one example for which such a radius measurement 
provides the identification of its nature. This is the case 
of Hat-P-2b, a 8 Mj mass ob ject, closely orbiting an F 
type star fsee lWinn et al.ll2007l and Table [TJ). As illustrated 
in Fig. [51 an irradiated brown dwarf model (long-dashed 
line) overestimates the radius by ~ 5%. Models including a 
340 M ffi core masfl can explain the measured radius (dash- 
dotted line) . Although this amount of heavy material is 
about the limit of what is available for planet formation, ac- 
cording to current core-accretion models (a s estimated from 
Eq. (1 ) for the HAT-P-2 system; see also iMordasini et alj 
(2009)), the presence of such a metal enrichment (Z < 15%) 
provides the simplest plausible explanation for the observed 
radius of HAT-P-2b, according to the present theory. As 
mentioned earlier, this 340 core for HAT-P-2b should 
be seen as a rough estimate of the upper limit for the avail- 
able heavy material in the system, but this analysis shows 
that the currently observed low radius of this object can- 
not be explained without a substantial enrichment. It would 
certainly be interesting to see whether planet models from 
other groups yield or not similar determinations. Given the 
fact that these various planet models share many common 
physics inputs (in particular the H/He and heavy element 
EOS), it would be surprising that they reach severely dif- 
ferent conclusions. While keeping in mind the remaining 
uncertainties in planet cooling theory, the present analysis 
provides - with the parameters observed so far - a confirma- 
tion of the validity of the core-accretion model, and makes 
Hat-P-2b the first confirmed 8 Mj genuine planet formed 
by core-accretion in a proto-planetary disk. 

6. Conclusion 

This work focusses on the modelling of the first confirmed 
transiting planets discovered by the CoRoT mission. We 
have distinguished three sorts of objects. First, planets 
whose radius can be explained by standard structure and 
evolution planetary models including the effect of irradia- 
tion from the parent star. For these objects, such as CoRoT- 
Exo-4b, there is no need to invoke extra physical mecha- 
nisms and, in that case, an upper limit can be inferred on 
their global heavy material content. The second category 
of objects is characterised by an abnormally large radius, 
with two examples in the current CoRoT sample. Ground- 
based surveys already found a significant fraction of such 

1 Note that - as discussed in jj2]- this amount of heavy material 
does not necessarily need to be in a core but could be distributed 
all over the planet. 
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Fig. 5. Hat-P-2b. Long dashed line: Cooling sequence 
of an irradiated 8.04Mj brown dwarf. Dash-dotted line: 
Irradiated case with a 340 M e core. Red box: observational 
la error bar. 



planets and CoRoT confirms this trend. We show that a 
small but finite eccentricity of ~ 0.02 provides enough tidal 
energy dissipation to explain the radius of CoRoT-Exo-lb. 
For CoRoT-Exo-2b, a significantly larger value, e ~ 0.15, 
is required. We emphasize the fact that a zero eccentricity 
value is assumed in the light curve and radial velocity data 
analysis of current CoRoT planets, and of many other tran- 
siting objects. This hypothesis is based on the idea that (i) 
tidal circularisation of the orbit is the asymptotic equilib- 
rium state of these planets, (ii) this circularisation occurs 
on short timescales compared to the system's age. Two as- 
sumptions wh ich have recent l y been shown not to be neces- 
sarily correct (jLevrard et a l. (2009)). Eccentricity determi- 
nations of transiting planets, when unknown, should thus 
be redone in light of these results. Our estimates for the 
required eccentricity to lead to enough tidal dissipation to 
explain the observed radii could thus be compared with 
future follow-up observations/determinations. 

The third category of planets includes the "massive" 
object CoRoT-Exo-3b, in the overlapping mass regime be- 
tween brown dwarfs and giant planets. We show that the 
remaining large uncertainties in the radius determination 
(~ 7%) for this object do not allow a clear identification of 
its nature. For this object's mass, present models predict 
less than 3% difference between the radius of a brown dwarf 
(solar composition) and of a planet with a realistic heavy 
material enrichment. By itself, the agreement between the 
observed radius and the theoretical predictions is a beauti- 
ful confirmation of the validity of these latter in the brown 
dwarf mass range, and brings confidence in the validity of 
the physics included in these models. Although our analysis 
is inconclusive for CoRoT-Exo-3b, given the present radius 
uncertainties, we show that the radius determination is dis- 
criminant in the case of Hat-P-2b, which is thus the first 
confirmation of the possibility to form massive planets by 
core accretion (possibly with subsequent collisions) up to 
m>8Mj. Finally, our analysis shows that, according to 
the present models, a typical < 5 % accuracy on the radius 
determination must be achieved in future space-based or 
ground based transit detections to clearly distinguish plan- 
ets from brown dwarfs in their overlaping mass domain. We 
stress, of course, that the conclusions of the present paper 
are based on presently published observational error bars 
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and on our present theoretical description of planet struc- 
ture and evolution, still hampered by many uncertainties. 
The present results thus need to be confronted to e.g. future 
COROT or KEPLER detections in order to assess their re- 
liability. 
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